First observation of the hyper superheavy hydrogen ^H 
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Abstract 

Three candidate events of the neutron-rich hypernucleus ^H were uniquely identified in the 
FINUDA experiment at DA$NE, Frascati, by observing 7r+ mesons from the [K^^ ,it^) pro- 
duction reaction on ^Li targets, in coincidence with 7r~ mesons from ^H — > ^He-f tt^ weak decay. 
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Details of the experiment and the analysis of its data are reported, leading to an estimate of 
(2.9 ±2.0) • 10^^ / K^ for the ^H production rate times the two-body tv" weak decay branching 
ratio. The ^H binding energy with respect to ^H + A was determined jointly from production 
and decay to be Ba ~ (4.0 ± 1.1) MeV, assuming that ^H is unbound with respect to '^H + 2n 
by 1.7 MeV. The binding energy determined from production is higher, in each one of the three 
events, than that determined from decay, with a difference of (0.98 ± 0.74) MeV here assigned 
to the Og.s. — >■ 1^ excitation. The consequences of this assignment to A hypernuclear dynamics 
are briefly discussed. 
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1. Introduction 

The role of the A hyperon in stabilizing nuclear cores was pointed out long ago by Dalitz 
and Levi Setti [1] as part of a discussion focusing on light hypernuclei with large neutron 
excess. This property is demonstrated by the observation of ^Hc, Ijie, ^He, ^Be and 
^J(B hypernuclei in emulsion experiments [2]. No unstable-core hydrogen A hypernuclei 
have been established so far, although the existence of the lightest possible one ^H was 
predicted by Dalitz and Levi Setti [1] and subsequently reinforced in estimates by Majling 
[3]. The neutral-baryon excess in ^H, in particular, would be {N + Y)/Z ~ 5, with Y = 1 
for a A hyperon, larger than the maximal value in light nuclei, N/Z = 3 for *He [4]. 
Neutron-rich light hypernuclei could thus go beyond the neutron drip line for ordinary 
nuclear systems. 

Two-body reactions in which neutron rich hypernuclei could be produced are the fol- 
lowing double charge-exchange reactions: 

K-+^Z^i{Z-2)+7r+, (1) 

induced on nuclear targets by stopped K^ mesons or in flight, and 

n-+^Z^i{Z-2) + K+ (2) 

with 7r~ mesons in flight (p^- > 0.89 GeV/c). 

The simplest description of the above reactions is a two-step process on two different 
protons of the same nucleus, converting them into a neutron and a A, with the additional 
condition that the final nuclear system is bound. For (1) it amounts to K^p — > Att'' re- 
action followed by 7r"p -^ nn^ or K^p — > K'^n followed by K'^p -^ Att^, for (2) to a 
TT^p — > riTT*^ reaction followed by ir'^p — > K^ A or n^p — > K'^A followed by K'^p — > K^n. 
Another mechanism is a single-step double charge exchange m~p — >■ Tj~m't (where m 
stands for meson) feeding the S component coherently admixed into the final A hyper- 
nuclear state. Such admixtures are essentially equivalent to invoking a second step of 
S^p — > An conversion. These two-step processes are expected to occur at a much lower 
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rate (reduction factor < 10^^ [5]) than the production of normal A hypernuclei by means 
of the corresponding single-step two-body reactions {K^,tt^) and {tt^,K~^). 

The first experimental attempt to produce neutron-rich hypernuclei by the reaction 
(1) with K^ at rest was carried out at KEK [6]. Upper limits were obtained for the 
production of ^He, ^^Be and ^^C hypernuclei (on ^Be, ^^C and ^^O targets respectively) 
in the range of (0.6 — 2.0) • lO^'^/K^^ , while the theoretical predictions for ^^Be and 
^^C [7] lie in the interval (10^^ — 10^''^) /K^, which is at least one order of magnitude 
lower than the experimental upper limits and three orders of magnitude smaller than the 
standard one-step {K~^^p,TT~) reaction rates on the same targets {10^^ / K~^.^p) . 

Another KEK experiment [8] reported the observation of ^^Li in the {tt^ ,K^) reaction 
on a ^°B target with a 1.2 GeV/c 7r~ beam. A production cross section of 11.3 ± 1.9 
nb/sr was evaluated; the result, however, is not directly comparable with theoretical 
calculations [9] since no discrete structure was observed and the production cross section 
was integrated over the whole bound region (0 < J5a < 20 MeV). 

A further attempt to observe neutron-rich hypernuclei by means of the reaction (1), 
with K^ at rest, was made at the DA$NE collider at LNF by the FINUDA experiment 
[10], on ^Li and ''Li targets. The limited data sample collected during the first run 
period of the experiment was used to estimate the production rates per stopped K^ 
of ^H and ^H. The inclusive tt"*" spectra from ^Li and ^Li targets were analyzed in 
momentum regions corresponding, through momentum and energy conservation, to B\ 
values discussed in the literature. Because of the dominant contribution of the reactions 



^stop + P ^ S+ + 7r" 



n 



+ TT+ (- 130 < p^+ < 250 McV/c) (3) 



and 

^n + 7T+ (- 100 < p^+ < 320 McV/c), (4) 

which give the main component of the inclusive 7r+ spectra for absorption of stopped K'^ 
mesons on nuclei, and owing to a limited statistics, only upper limits could be evaluated 
for A hypernuclear production: 

i?^+ (iH) < (2.5 ± 0.4,tatiE!1syst) • 10- V^stop, (5) 

K^ (iH) < (4.5 ± 0.9,tat+l];tsy.t) • 10- V^stop, (6) 



in addition to an upper limit determined in ^^C: 

K+ ClBe) < (2.0 ± 0.4,tattl!:Lyst) • 10- Vi^stop, (7) 

lowering by a factor ^ 3 the previous KEK determination [6] . 

In this article we present the analysis of the total data sample of the FINUDA exper- 
iment, collected from 2003 to 2007 and corresponding to a total integrated luminosity 
of 1156 pb^^, aiming at assessing the existence of ^H and determining the production 
rate by means of the {K^^ jir^) reaction on ®Li targets. A preliminary account of the 
results, reporting three clear events of ^H, appeared in [11]. 
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Fig. 1. Left: Binding energy scheme for a system of one proton, four neutrons and one A relative to the 
summed mass of ^H + A, M = 5805.44 MeV, with the (blue in the web version) hatched box denoting 
the width of ^H [12]. The two lowest horizontal lines stand for predictions from Refs. [1,13]. Right: mean 
value of the ^H g.s. mass obtained in the present analysis jointly from production and decay; the (red 
in the web version) shaded box represents the error on the mass mean value obtained from the three ^H 
events reported here. 

The binding energy of ^H with respect to the unstable ^H core was estimated in 
Refs. [1,3] as B\ = 4.2 JVIeV, making ^H particle stable with respect to its ^H + 27i lowest 
threshold, as shown in Fig. 1. We recall that the binding energy Ba of hypernucleus ;^Z 
is defined as: 

Ba = Afcoro + Ma -Ma^, (8) 

where Mcorc is the mass of the ('^^^^Z core nucleus in its ground state (g.s.), as deduced 
from the atomic mass tables [14]. The ^H nuclear core, colloquially termed "superheavy 
hydrogen", was observed as a broad resonance (1.9 MeV FWHM) at energy about 1.7 
MeV above the '^H + 2n threshold [12]. A substantially stronger binding, B^ijji) — 5.8 
MeV, was predicted by Akaishi et al. [13] for the 0+ g.s. on the basis of a coherent AiV — 
Y,N mixing model originally practised for the ^H cluster [15]. This coherent AA^ — Y.N 
mixing induces a spin-dependent ANN three-body interaction which affects primarily 
the 0+ g.s., increasing thus the « 1 MeV 1+ excitation expected from ^H to 2.4 MeV in 
^H. If this prediction is respected by Nature, it could imply far-reaching consequences 
to strange dense stellar matter. 

In the next sections we describe briefly the FINUDA experimental apparatus, and the 
analysis technique applied to the data collected on ^Li targets. We then report on three 
^H candidate events found by observing 7r+ mesons from production and tt^ mesons 
from decay in coincidence. These events prove robust against varying the cuts selected 
in the analysis, and give evidence for a particle stable ^H. The measurement background 
is evaluated and the production rate of ^H is estimated. We end with a brief discussion 
of the ^H excitation spectrum as constrained by the three candidate events. 



2. Experimental apparatus 

FINUDA was a hypernuclear physics experiment installed at one of the two interac- 
tion regions of the DA$NE e+e" collider, the INFN-LNF $(1020)-factory. A detailed 
description of the experimental apparatus can be found in Ref. [16]. The layout figured 
a cylindrical symmetry arrangement; here we briefly sketch its main components moving 
outwards from the beam axis: the interaction/target region, composed by a barrel of 12 
thin scintillator slabs (TOFINO), surrounded by an octagonal array of Si microstrips 
(ISIM) facing eight target tiles; the tracking device, consisting of four layers of position 
sensitive detectors (a decagonal array of Si microstrips (OSIM), two octagonal layers of 
low mass drift chambers (LMDC) and a stereo system of straw tubes (ST)) arranged in 
coaxial geometry; the external time of flight detector (TOFONE), a barrel of 72 scintil- 
lator slabs. The whole apparatus was placed inside a uniform 1.0 T solenoidal magnetic 
field; the tracking volume was immersed in He atmosphere to minimize the multiple 
scattering effect. 

The main features of the apparatus were the thinness of the target materials needed 
to stop the low energy (~ 16 MeV) K~^s from the $ -^ K~K'^ decay channel, the 
high transparency of the FINUDA tracker and the very large solid angle (~ 27r sr) 
covered by the detector ensemble; accordingly, the FINUDA apparatus was suitable to 
study simultaneously the formation and the decay of A hypernuclei by means of high 
resolution magnetic spectroscopy of the emitted charged particles. 

In particular, for 7r+ with momentum ~ 250 MeV/c the resolution of the tracker can be 
evaluated by measuring the width of the momentum distribution of the monochromatic 
(235.6 MeV/c) /i"*" coming from the Kfj_2 decay channel; for reactions occurring in the 
apparatus sector where ^Li targets were located, it is a-p = (1.1 ± 0.1) MeV/c [17]; the 
precision on the absolute momentum calibration, obtained from the mean value of the 
same distribution, is better than 0.12 MeV/c for the ^Li targets, which corresponds to a 
maximum systematic uncertainty in the kinetic energy axsysiT^^) =0.1 MeV. 
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Fig. 2. Distribution of low momentum tt~ from ®Li targets. The continuous black curve represents a fit 
the the spectrum given by the sum of a second degree polynomial [dashed (red in the web version) curve] 
and a gaussian function [dot-dashed (blue in the web version) curve]. For more details, see text. 



For TT^ with momcntuni ^ 130 MeV/c the resolution and absolute calibration can be 
evaluated from the momentum distribution of the monochromatic tt^ coming from the 
two-body mesonic weak decay of ^^H, produced as hyperfragment with a formation prob- 
ability of the order of 10^^ — 10^^ per stopped K^ [18]. Figure 2 shows the distribution 
for low momentum tt^ from ®Li targets, before acceptance correction; the spectrum is 
fitted in the 120-140 MeV/c momentum range (continuous black curve) with the sum of 
a second degree polynomial, representing the background from quasi-free A decay and 
quasi- free E'^ production (dashed (red in the web version) curve in the figure), and a 
gaussian function representing the ^H mesonic decay contribution (dot-dashed (blue in 
the web version) curve); the fit gives a x^/ndf = 79.1/74, a mean Hp = (132.6 ± 0.1) 
MeV/c and a standard deviation dp — (1.2 ± 0.1) MeV/c for the gaussian function, 
directly measuring the experimental resolution. For comparison, p^- = (132.80 ± 0.08) 
MeV/c from i?A(^H) — 2.04 ±0.04 MeV, as determined from emulsion studies [2]; hence 
the absolute uncertainty is 0.2 MeV/c. and the corresponding systematic uncertainty in 
the kinetic energy is then crTsj/s(7r~) = 0.14 MeV. 

To perform particle identification, the information of the specific energy loss in both 
OSIM and the LMDC's is used; the mass identification from the time of flight system 
(TOFINO-TOFONE) for high momentum tracks is also used. The final selection is per- 
formed by requiring the same identification from at least two different detectors. 

3. Analysis technique 

In the second data taking the statistics collected with ^Li targets was improved by 
a factor 5 with respect to the first run. However, even with the improved statistics, we 
could not observe in the inclusive 7r+ spectra clear peaks that could be attributed to the 
two-body reaction: 

K-,^^ + ^U^lli + 7r+ (p,+ ^ 252 MeV/c). (9) 

Exploiting the increased statistics, we tried then to reduce the background overwhelming 
the events from reaction (9) by examining the spectra of 7r+ in coincidence with the tt^ 
coming from the mesonic decay of ^H: 

^H -^ ^He + TT- {p^- - 130 - 140 MeV/c). (10) 

The branching ratio for (10) is expected to be about 50% taking into account the value 
measured for the analogous decay ^H — > ^He -I- tt" [18]. (7r+, tt") coincidence events, 
associated with K^^s stopped in the ^Li targets, were thus considered; only reaction (3) 
contributes to the background of this sample. 

We examined thus the two-dimensional raw spectrum of 7r+ versus n^ momentum, 
shown in Fig. 3, in order to recognize possible enhancements due to occurrence of the re- 
actions (9) and (10) in sequence. The low statistics and the strong background prevented 
us from finding statistically significant accumulations of events in the plot arising from 
a bound ^H. 

In order to isolate the events due to the possible formation of a bound ^H, we con- 
sidered energy conservation for both reactions (9) and (10). Momentum conservation is 
automatically ensured by the fact that both reactions of formation (9) and decay (10) 
occur at rest. The stopping time of ^H in the material is indeed shorter than its lifetime. 
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Fig. 3. 7r+ momentum vs tt inomentuin for ''Li targets. The (red in the web version) band dots stand 
for events with T{tt+) + T(tt^) = 202 - 204 MeV. See text for more details. 

For (9) we may write explicitly: 

M{K-) + UI{p) + 3M{n) - B^U) = Af (^H) + THE) + M{n+) + T{tt+), (11) 

in which, in obvious notation, M stands for a particle mass, T - its kinetic energy, and 
S(^Li) - the binding energy of ®Li. For (10) we may write: 



M(^H) = 2M{p) + 4M(n) - B(^Hc) + T{^nc) + Mi-R-) + T{tt-), 



(12) 



in the same notation as above. Combining Eqs. (11) and (12) in order to eliminate 
M (^H), we get the following equation: 



T(7r+)+T(7r~) 



: M{K^) + M{p) - M{n) - 2M{tt) 
-B(^Li) + B{^Yie) - T(^He) - r(^H). 



(13) 



All the terms on the right-hand side are either known constants or quantities that can 
be evaluated from momentum and energy conservation, except for T(^H) (T(^He)) that 
depends explicitly (implicitly) on the unknown value of i?A(^H). A variation of iJA^H) 
between and 6 MeV introduces a change of ^ 0.3 MeV in the kinetic energy T(^H) 
in (11), corresponding to a sensitivity of 50 keV per MeV of Bf^Q^)^ and a change of 
~ 0.2 MeV in T(7r+) + T{'k~) in (13), corresponding to a sensitivity of 30 keV per MeV 
of Bfs^lj^). These variations are much lower than the experimental energy resolutions 
for 7r+(250 MeV/c) and 7r-(130 MeV/c): crT(7r+) = 0.96 MeV and ariTT-) = 0.84 MeV. 
The FINUDA energy resolution for a (tt+jTt^) pair in coincidence is therefore ctt = 
1.3 MeV, where cttcxp = ^0.96^ + 0.842 ^ ^ 3 j^eV is the total 

0.17 MeV is 
5 



Texp "T "Tsys 

experimental energy resolution and o'T'sys 



V^O-Tsys ( 7r+ ) 2 + crygy s ( TT ) 2 

the total systematic error on energy. To be definite, we assume a value of i3A(^H) 
MeV, halfway between the conservative estimate of 4.2 MeV [1,3] and Akaishi's prediction 
of 5.8 MeV [13]. The r.h.s. of Eq. (13) assumes then a value of r(7r+)+T(7r-) = 203.0±1.3 
MeV. 
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Fig. 4. Distribution of raw total kinetic energy T(7r+) + T{n ) for (7r+, n ) coincidence events from ^Li 
targets. The (red in the web version) shaded vertical bar represents the cut T(7r+) + T{it~) = 202 — 204 
MeV. 



We considered then the raw spectrum of the total kinetic energy, T{tt^) + T(7r~), for 
the coincidence events, shown in Fig. 4. Events in the summed energy distribution were 
selected in the region (203 ± 1) MeV, indicated by the (red in the web version) filling 
in the figure. The half-width of the interval corresponds to ~ 77% of the FINUDA total 
energy resolution; the value was chosen as a compromise between the strong requirement 
of reducing the contamination from background reactions, as will be discussed in more 
detail in the following, and the plight for reasonable statistics, leading to application of a 
selection narrower than the experimental resolution. The selected events are represented 
by red dots in Fig. 3. 

The raw distributions of p^+ and p^^- for the events selected are shown in Fig. 5 by 
the continuous line histogram, falling off to zero at p^+ = 245 MeV/c in the higher 
momentum region, and at p^^- = 145 MeV/c in the lower momentum region. These 
limiting values, when inserted in Eqs. (9) and (10) for two-body kinematics ^H production 
from rest and decay at rest, yield ^H mass values higher than the total mass of both 
(A + ■^H + 2n) and (A + ^H) thresholds marked in Fig. 1. A ^H mass equal to the mass 
of its lowest particle stability threshold ^H + 2n corresponds to values of p^+ = 251.9 
McV/c and p„- — 135.6 MeV/c. A genuinely bound ^H system, therefore, requires that 
pion momenta satisfying p^+ > 251.9 MeV/c and p^- < 135.6 MeV/c are selected. 
The cuts actually applied in the analysis of the data, p^+ = (250 — 255) MeV/c and 
Pjt- = (130 — 137) MeV/c, as marked by the (blue in the web version) shaded vertical 
bars in Fig. 5, allow for a wide range of ^H masses from the (A + ''H + 2n) threshold, 
about 2 MeV in the ^H continuum, down to B\{\}i) > 6 MeV, somewhat below the mass 
predicted by Akaishi [13]. These cuts do not exclude completely an eventual contribution 
from the production and decay of ^^^H + 2n, of a weight which is anyway negligible, as 
discussed in the next section. 
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Fig. 5. Continuous histograms: distribution of 7r+ (upper part) and tt" (lower part) momenta for (7r+ , n^ ) 
coincidence events with T(7r+) + T(7r~) = 202 — 204 MeV from ®Li targets before acceptance correction. 
The (blue in the web version) shaded vertical bars indicate the final selection regions. Dashed (red in 
the web version) histograms represent the S+ background spectra: see text for more details. 

4. Results 



Three events, out of a total number of ~ 2.7-10^ K detected at stop on the ^Li targets, 
satisfy the final requirements, r(7r+) +T(7r-) = 202-204 MeV, p^+ = 250-255 MeV/c 
and Pt^- — 130 — 137 MeV/c. These events, within the (red in the web version) shaded 
rectangle on the l.h.s. of Fig. 6, are candidates for ^H. The 7r+ momenta which this 
rectangle encompasses go up from a value corresponding to the (A + "^H + 2ti) threshold 
to a value corresponding to the binding energy predicted by Akaishi, whereas the tt" 
momenta which the rectangle encompasses go down from a value corresponding to the 
same (A + ■^H + 2n) threshold to about 2cr(p^-) below the value predicted by Akaishi 
[13]. 

Different choices of T(7r+) + T{'k ) interval widths (2 — 6 MeV) and position (center 
in 202 - 204 MeV) and of p^+/p^- interval widths (5 - 10 and 8-15 MeV/c) with 
fixed limits at 250 and 137 MeV/c respectively to exclude the unbound region, affect the 
populations of the corresponding single spectra but not the coincidence spectrum. As 
an example, in Fig. 6 a comparison is made between the (Ptt+jPtt-) plots satisfying the 
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Fig. 6. tt"*" momentum vs n momentum for ^Li target events with T{n~^) + T{it ) = 202 — 204 MeV 
(l.h.s.) and with T{tt+)+T{tt^) = 200-206 MeV (r.h.s.). The shaded (red in the web version) rectangles 
on each side consist of a subset of events with p^+ = 250 — 255 MeV/c and p^- = 130 — 137 MeV/c. 
The hatched (blue in the web version) rectangles on each side are symmetric subsets of events to those 
in the shaded rectangles. 



actual selection T(7r+) + T{Tr^) = 202 - 204 MeV (l.h.s.), and similar plots admitting a 
wider selection range T(7r+) + T{-k^) = 200 — 206 MeV (r.h.s.). The global population 
increases for the wider cut, as expected, but the events that satisfy simultaneously also 
the separate selections imposed on p„+ and p^r- (shaded rectangles in the upper left part 
of the plots) remain the same. A similar stability is not observed in the opposite corner 
of the plots where, on top of the events already there on the left plot, five additional 
events appear on the right plot upon extending the cut. Quantitatively, fitting the pro- 
jected TT* distributions of the l.h.s. of Fig. 6 by gaussians, an excess of three events in 
both p^± distributions is invariably found, corresponding to the shaded (red in the web 
version) rectangle. The probability for the three events to belong to the fitted gaussian 
(background) distribution is less than 0.5% in both cases. It is possible, moreover, to see 
directly from the two dimesional plots that variations of the independent momentum se- 
lections do not produce any effect. Systematic errors due to the applied analysis selection 
are thus ruled out. 

It is also worth noticing that the tight momentum cuts imposed on the (7r+,7r^) coin- 
cidence events allow to eliminate completely any contamination due to possible 7r^/e^ 
misidentification. Furthermore, /i^'^'s from iiTJ'j decay are clearly separated from tt'^'s 
coming from the opposite K^ interaction vertex. Figure 7 shows a front view of one of 
the three events, as reconstructed by FINUDA. 

By evaluating event by event the corresponding ^H mass from both production (9) 
and decay (10) reactions, the mass values listed in Table 1 are obtained. A mean value 
for each event jointly from production and decay was also evaluated, with an error given 
by the highest of the two uncertainties, 0.96 MeV. For the global mean mass value of 
the three events we then find A/(^II) — (5801.43 ± 0.55) MeV, where the uncertainty 
0.55 = 0.96/-\/3 MeV refiects the uncertainty assigned to each event. This global mass 
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Fig. 7. Left: front view of one of the ^H candidate events reconstructed by FINUDA where a (tt"*", n~) 
pair emerges from a ^Li target and crosses the spectrometer. Right: expanded view of the target region 
for the same event where the K~ track stops in a ®Li target. 

Table 1 

Kinematical properties and ^H mass, Af(^H), of the three ^H candidate events from production (9) and 
decay (10) reactions. Listed in the last two columns are the mean and the difference of the production 
and decay masses. Ttot indicates the total n^ kinetic energy T(7r+) + T{tt~). 



Ttot 
(MeV) 



202.6±1.3 
202.7±1.3 
202.1±1.3 



(MeV/c) 



251.3±1.1 
250.1±1.1 
253.8±1.1 



Pn- 
(MeV/c) 



135.1±1.2 
136.9±1.2 
131.2±1.2 



Af(iH) 
prod. (MeV) 



5802.33±0.96 
5803.45±0.96 
5799.97±0.96 



M(iH) 
decay (MeV) 



5801.41±0.84 
5802.73±0.84 
5798.66±0.84 



M(6h) AM(6H) 
mean (MeV) (MeV) 



5801.87±0.96 0.92±1.28 
5803.09±0.96 0.72±1.28 
5799.32±0.96 1.31±1.28 



uncertainty, however, is considerably smaller than the mean-mass spread of the three 
events. We therefore decided to relax the assigned uncertainty by calculating it from the 
spread of the three mean mass values, which yields 



Af(^H) = 5801.4 ± 1.1 MeV, 



(14) 



with uncertainty larger than the 0.96 MeV and 0.84 MeV measurement uncertainties in 
production and decay respectively. The standard deviation of this uncertainty is 0.55 
McV which together with cr = 1.1 MeV is still short of the 2.11 MeV deviation of the 
third event mean mass from the mean mass value. This observation could indicate some 
irregularity in the reconstruction of the third event. To regain confidence, each one of the 
three events was checked visually for irregularities but none was found. The third event, 
in particular, is shown in Fig. 7. 

Listed in the last column of Table 1 are values of AA/(^H), defined as the difference 
between the ^H mass values obtained from production and from decay. The mass values 
obtained from production are systematically higher than those from decay by 
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AM(^H) = 0.98 ± 0.74 MeV, (15) 

where the uncertainty is evaluated from the 1.3 McV uncertainty for T(7r+) + T{tt^) 
from which each of the mass differences is directly determined. Unlike the mean M(^H) 
mass value, the spread of the production vs decay mass differences is well within la. A 
possible physical origin of the AM(^H) systematics is discussed in a subsequent section. 
The mean mass value (14) corresponds to a ^H binding energy Ba — (4.0 ± 1.1) MeV 
with respect to the (A + ^H) threshold, and to Ba = (0.3 ± 1.1) MeV with respect to the 
lowest threshold (^H + 2n). The ^H mean mass value and its uncertainty are indicated 
on the r.h.s. of Fig. 1 with respect to the various thresholds and predictions shown on 
the l.h.s. of the figure. 

5. Background estimation 

Before discussing the physical interpretation of the above results, it is mandatory to 
check carefully that the three observed events do not arise from physical or instrumental 
backgrounds that could affect the data. Concerning the physical backgrounds, a complete 
simulation has been performed of possible ^j^op absorption reactions on both single 
nucleons and pairs of strongly correlated nucleons that lead to the formation and decay 
of A and S hyperons. Of these reactions, only the following chain leads to (tt+jTt") 
coincidences in the same momentum ranges corresponding to the production and mesonic 
decay of ^H and which are respected by the three candidate events: 

X.top + ®Li ^ S+ + ''He + n + TT- {p^- < 190 MeV/c) 

^n + Tr+ {p^+ < 282 MeV/c). (16) 

This reaction chain has been studied by means of the FINUDA simulation program 
fully reproducing the apparatus geometry, detection efficiency and the trigger efficiency. 
The interaction of K~ with the target nucleus has been simulated with two different 
approaches. In the first approach, the quasi-free approximation was adopted for the 
interaction of the K^ with a proton of the target nucleus, K^^. + p ^ S+ + tt^, 
taking into account the nucleon Fermi motion; the residual nucleus was considered as a 
spectator and the notation "^He + n" is just a label to indicate that the system is highly 
particle unstable. Pions arising from (16) were processed by the pattern recognition and 
reconstruction programs of FINUDA as real data. In common with all simulated reaction 
chains, the simulated events were then submitted to the same quality cuts and to the 
same selections criteria applied in the data analysis. Three events were found out of a 
total of 2.2 • 10^ K^ mesons simulated to stop on ^Li targets and forced to undergo the 
(16) "quasi- free" reaction chain with a probability of 1. Taking into account the number 
of actual K^ mesons stopped on ^Li targets, the branching fraction for the ^i^op +P ^^ 
S+-|-7r^ reaction on nuclei measured on ^^C [19] and on ''He [20], (0.159±0.012) evaluated 
as a weighted mean, the S+ + n ^ A + p conversion probability [21], (0.45 ± 0.04), and 
the S+ ^ n + 7r+ decay branching ratio, (0.483 ± 0.003), an expected E+ background of 
0.15 ± 0.09 events on ^Li targets is obtained. 

In a second approach, the interaction of K'^ mesons with the target nucleus as a whole 
was considered, applying directly the 4-body kinematics to (16). Five events were found 
out of a total of 2.7T0^ K~ mesons simulated to stop on ^Li targets and forced to undergo 
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the (16) "4-body" reaction chain with a probabiUty of 1. Taking into account the same 
normahzation factors used for the "quasi-free" approach, an expected S"*" reaction chain 
background of 0.20 ± 0.11 events on ^Li targets was obtained under the hypothesis that 
S+ production on ^Li in this approach always gives a recoihng ''He nucleus. Final states 
corresponding to further fragmentation of the ^Li target nucleus, such as K~^ + ^Li — >■ 
E+ + '^He + n + n + tt~ , give weaker background contribution, owing to the requirements 
imposed on T(7r+) + T{tt^) (< 180 MeV for final states of the E+ production reaction 
with more than 4 bodies) and on the 7r+ and tt^ momenta. 

We also considered the distortion of Eq. (16) reaction chain spectra due to the '^He + n 
final state interaction leading to ^He, a resonance centered at ^ 0.8 MeV above the 
''He + n threshold with F = 1.36 MeV [22]. To this end we required that once the ''He 
and neutron momenta generated by the 4-body phase space simulation corresponded to 
the formation of the ^He resonance, the momenta of the remaining particles, E"*" and 7r~, 
should be modified accordingly. We passed then these modified phase space distributions 
through the selection criteria described above and found variation of less than 1% in the 
background value evaluated for a 100% ''He + n final state. 
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Fig. 8. T(7r+) + T{iT~) distributions; (red in the web version) points: experimental data; (blue in the 
web version) dashed histogram: "quasi-free" simulation; (violet in the web version) dotted histogram: 
"4-body" simulation; (black) continuous histogram: best fit to the data with fractions of the simulated 
templates. The simulated distributions have been normalized to the experimental distribution area. For 
more details see text. 

In Fig. 8 the experimental r(7r+)-|-r(7r^) spectrum is shown together with spectra ob- 
tained from the "quasi-free" and "4-body" simulations of the (16) process: the simulated 
spectra were normalized to the area of the experimental distribution. As may be seen, 
the "quasi-free" spectrum (dashed (blue in the web version) histogram) reproduces the 
experimental distribution better than the "4-body" (dotted (violet in the web version) 
histogram), but exhibits a too sharp decrease in the 200-210 MeV region and underesti- 
mates the low energy tail. To obtain a satisfactory description, a fit of the experimental 
spectrum was performed with fractions of the two simulated templates; a standard likeli- 
hood fitting method, using Poisson statistics, was applied in which both data and Monte 
Carlo statistical uncertainties were taken into account [23]. Particular care was devoted 
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to the description of the 200-210 McV slope, where the selection of the ^H candidate 
events is made. The continuous (black in the web version) histogram in Fig. 8 represents 
the best fit to the 180-220 MeV region; the resulting fractions are 0.743 ± 0.019 and 
0.257 ± 0.017 for the "quasi-free" and "4-body" templates respectively, with a x^/ndf = 
40.0/39. We note that varying the width of the fit region from 180-220 MeV to 130-220 
MeV spoils the fit, increasing the x^/ndf value by a factor of ^ 3.8, while the fractions 
of the two templates change by less than 0.025, corresponding to 1.3 — 1.5 a. 

Back in Fig. 5, the dashed (red) histograms represent the separate p^+ and p^- dis- 
tributions obtained by adopting the above fractions to the events successfully simulated 
within the two approaches discussed above and satisfying the cut T{tt^) + T{tt^) = 
202 — 204 MeV. Only a qualitative agreement is reached with the experimental, low 
statistics distributions. The estimated background spectra are shifted from the experi- 
mental ones toward higher momentum for the n^ spectrum and toward lower momentum 
for the 7r~ spectrum. The difference is significant over the statistical fiuctuation. How- 
ever, these shifts cause overestimates of the background counts when the tails of the 
simulated distributions are used for an estimate. It is thus possible to conclude that the 
background estimate is safe in spite of the slight deviations noted here. In particular, in 
the ^H selected regions, indicated by shaded (blue in the web version) vertical bars, the 
contribution due to the S^ background for events satisfying also the cuts on p^+ and 
Ptt- in coincidence corresponds to 0.16 ± 0.07 events on ^Li targets (BGDl). 

Another reaction chain capable of providing background events is 

X.top + ®Li ^ ^H + 2n + 7r+ (p„+ < 252 McV/c) 

Hc + TT- (p„- - 132.8 MeV/c). (17) 



stop 



The momentum of the ^H decay tt^ is close to the momentum of the tt^ from the two- 
body decay of ^H, p^- ~ 134 MeV/c, evaluated assuming B\ = 5 McV which is halfway 
between the two theoretical estimates exhibited in Fig. 1. 

The probability of having background contribution from this reaction chain was eval- 
uated taking into account the phase space fraction of the reaction (17) available for tt+'s 
satisfying the momentum selection p^+ = 250 — 255 MeV/c, 4 • 10^^, and the probabil- 
ity for a K~f. to produce a ^H accompanied by a 7r+ on a ®Li target. In Ref. [18] the 
probability of producing ^H on ^Li targets was reported to be (3.0 ± 0.4) • 10^^ / K^^. 
and, furthermore, the probability of producing it together with a charged pion was in- 
dicated to be 0.49 ± 0.08. Using these values, the formation probability of ^H + tt* on 
®Li target, the closest isotope of ^Li, was evaluated to be (1.47 ± 0.15) • 10^'^ /K^^.. In 
addition, a branching ratio 0.49 [18] for the two-body decay ^H -^ ^He -I- tt^ has to be 
included. A total probability of (2.8 ±0.5) • 10^^ is obtained. From this value, taking into 
account the global efficiency of FINUDA (acceptance, reconstruction and analysis cuts) 
it is possible to evaluate a background of 0.05 ±0.01 reconstructed events for the 2.7- 10^ 
K~ mesons stopped on the ^Li targets. It should be noted that this value overestimates 
by far the actual contribution from (17) since the analysis of Ref. [18] is incapable of 
separating the contribution of ^H + 7r+ out of the global ^H + tt^ fraction. Note that 
a microscopic reaction evaluation of (17) requires a theoretical model which takes into 
account different channels (A production, E production, compound nucleus formation) 
the weights of which are not experimentally known. We chose to avoid relying on a model 
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and to assume, instead, an overly conservative evaluation of the background from the 
chain (17). The estimated level of this background is negligible with respect to the BGDl 
background from the chain (16). 

Other reaction chains have been considered, such as: 



(18) 



^stop + 'Li- 


-> S+ + iH + d + TT- {p^- < 165 McV/c) 




-^ n + 7r+ {p^+ < 250 MeV/c) 


i^stop + 'Li- 


-> in + 3n + 7r+ {p^+ < 242 MeV/c) 




^ ^Rc + n- {p^- - 115 McV/c) 


^stop + 'Li- 


-> A + ■''H + 2n + 7r+ {p^+ < 247 MeV/c) 



(19) 



^p + TT- {p„- < 195 McV/c) (20) 

These reaction chains may be safely discarded since all of them involve too low values of 
r(7r+) + T(7r~), less than 190 MeV, which are way outside the cut applied on Ttot(7i'); 
the chain (19) may be discarded in addition by the cut imposed on p^- . 

Finally, another mechanism which could produce a (7r+,7r^) pair in the final state is 

K-^^p + ^Li -^ ^Hc + 7r° {p^o - 280 MeV/c) 

^ *^Li + n- [p^- - 108 McV/c), (21) 

followed by a reaction on another ®Li nucleus: 

7r° + ^Li -^ ^He + 7r+ (p„+ - 280 MeV/c in the forward direction). (22) 

However, the kinematics of (21) rules out a contribution from this reaction chain when 
applying the cuts on tt^ momenta. Moreover, the mean free path of a i:^ with momentum 
~ 280 MeV/c, l^cc < 0.1 /xm, strongly reduces the probability of the second reaction of 
the chain with respect to 7r° decay. 

The main source of instrumental background could be the presence of fake tracks, due 
to fake signals from the detectors, misidentified as true events by the track reconstruc- 
tion algorithms. For this purpose we analyzed the events with a 7r+ and a 7r~ emitted in 
coincidence with a K~^ stopping in a given nuclear target under the following criteria: 

(i) events relative to target nuclei other than ^Li (^Li, ^Be, ^"^C, ^^O) were selected 
with the same selection criteria T{tt+) + T{tt-) = 202 - 204 MeV, p^+ = 250 - 255 
MeV/c and p^- = 130 — 137 MeV/c, as for ®Li. Incidentally, for events coming from 
the (16) reaction chain on nuclear targets heavier than ®Li, T(7r+) +T{Tr^) < 202 
MeV by at least 2(Tt so that this criterion actually selects the instrumental back- 
ground exclusively. Only one event was found, coming from ^Be target; 

(ii) events relative to the ®Li targets were selected with values of T(7r+) + T{tt^) — 
193 — 199 MeV so as to search for neutron-rich hypernuclei produced on the other 
targets. No events were found. 
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Taking into account the number of K^^ detected in ®Li targets and in all the other 
targets, we concluded that 0.27 ± 0.27 fake events should be expected from ®Li due to 
the instrumental background (BGD2). 

Combining together the expected number of events arising from physical and instru- 
mental backgrounds that affect our selected data, 0.16 ± 0.07 (BGDl) and 0.27 ± 0.27 
(BGD2), a total of 0.43 ± 0.28 background events has been established. From this value, 
following Poisson statistics, we may state that the three observed ^H candidate events do 
not belong to the background distribution with a confidence level of 99%; the difference 
between the measured yield and the total expected background can thus be safely consid- 
ered a ^H signal. The probability of observing three or more events from the background 
fluctuation following Poisson distribution with /i = 0.43 (BGD1+BGD2) or jj. — 0.16 
(only BGDl) is 0.0096 or 0.0006, respectively. In terms of a statistical significance S 
defined by 5 = C/VBGD, with a signal C = 3 — BGD, the statistical significance of the 
signal is 3.9 or 7.1, respectively. 

6. Production rate evaluation 

Using the background estimates of the last section, it is possible to evaluate the product 
R ■ BR(7r~), where R is the ^H production rate per stopped K~ and BR(7r~) is the 



branching ratio for the two-body weak decay ^H -^ ^Hc 



^ • B^^'^") - ^^^^f"! J'^m-^ = (1-3 ± 0-9) • lO^Vi^^op. (23) 

e(7r+) e{TT ) K^^^p{^U) 

In Eq. (23), e(7r+) and e(7r^) indicate the global efficiencies for 7r+ and tt^, respectively, 
including detection efficiency, geometrical and trigger acceptances and pattern recogni- 
tion, reconstruction and selection efficiencies, all of which have been evaluated by means 
of the full FINUDA simulation code, well tested in calculations for other reactions in 
similar momentum ranges [17,24,25]. K^ (^hi) is the number of K~ detected at stop in 
^Li targets. 

The value (23) has to be corrected for the purity of the ®Li targets used, 90%, for the 
0.77 <tt cut applied to T{tt^) + T(tt^), and for the fraction of ^H decaying in flight. 
In Ref. [18] a contribution of 20% is reported for the decay in flight of ^H produced on 
a ^Li target; extending this value also to ^H and considering that the cut applied to 
Ptt-, 130-137 MeV/c, allows to accept about one half of the pions emitted in flight, a 
correction factor of 10% is evaluated. The corrected result is 

[R ■ BR(^-)]eorr. = (2.9 ± 2.0) • 10-V^sTop- (24) 

By assuming BR(7r^) = 49%, in analogy to the weak decay ^H -^ *Hc + tt^ [18], 
we flnd R — (5.9 ± 4.0) • 10^^ / K^^. , fully consistent with the upper limit (5) obtained 
previously by FINUDA [10]. Although no theoretical calculation of this capture rate has 
been reported to provide direct comparison, the order of magnitude of the rate determined 
here is compatible with the interval of values calculated for production of heavier neutron- 
rich hypernuclei [7] with stopped K~ mesons and, as expected, is approximately three 
orders of magnitude lower than the capture rate for the production of ordinary particle- 
stable A hypernuclei. 
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7. Discussion 

The binding energy deduced from the three measured events listed in Table 1 was 
recorded in Eq. (14): -Ba(aH) = 4.0 ± 1.1 MeV with respect to ^H + A, close to the 
value BA(^He) = 4.18 ± 0.10 MeV with respect to ^He + A for the other known A = 6 
hypernucleus [2]. It is in good accord with the estimate 4.2 MeV made originally by 
Dalitz and Levi Setti [1] and confirmed by Majling [3]. It is lower by 1.8 MeV than 
the value 5.8 MeV suggested by Akaishi et al. [13], leaving little room for an attractive 
contribution from a ANN three-body force of a similar magnitude, 1.4 MeV, which in 
Akaishi's calculations arises from a coherent AiV — 'EN mixing model. This is consistent 
with a substantial weakening of AA^ — EN mixing contributions for the excess p shell 
neutrons in ^H with respect to the strong effect calculated in the s-shell hypernucleus 
^H [15]. Indeed, recent shell-model calculations by Millener indicate that AiV — EN 
mixing contributions to Ba and to doublet spin splittings in the p shell are rather small, 
about (10 ± 5)% of their contribution in ^H [26]. Nevertheless, given the measurement 
uncertainty of 1.1 MeV, one may not conclude that this ANN force contribution is 
negligible, but only that its influence appears considerably lower than predicted. For 
illustration, see Fig. 1. 

It is possible to avoid considering explicitly the AN— EN mixing effect in the evaluation 
of BA{%il) by updating the shell-model (SM) argument used in Ref. [1]. We adopt a 
cluster model for ^H in terms of ^H plus two p-shell neutrons coupled to J'^ = 0"^ as in 
^He g.s. The interaction of the A hyperon with this dineutron cluster, including any Ann 
force arising from A7V — EN mixing, may be deduced from A-He which consists of an a 
cluster plus precisely the same Ann configuration under consideration in ^H. Subtracting 
BA(iHe) = 3.12 ± 0.02 MeV from BaHRc), with a value -Ba(^Hc) = 5.36 ± 0.09 MeV 
obtained by extrapolating linearly from the known binding energies of the other members 
of the A = 7 hypernuclear T =1 isotriplet (see Fig. 3, Ref. [27]), we obtain 2.24 ± 0.09 
MeV for the Ann sum of two-body and three-body interactions involving the A hyperon. 
The value of BaHB.) is then obtained adding this 2.24 MeV to Ba{\E) = 2.04 ± 0.04 
MeV [2], so that B|^'(^H) = 4.28 ± 0.10 MeV.^ We have thus recovered the estimate 
originally made by Dalitz and Levi Setti [1] . 

As mentioned in the discussion of Table 1, the aH mass values obtained from pro- 
duction are systematically higher than the corresponding values obtained from decay, 
leading to a mass difference of AM(^H) = 0.98 ± 0.74 MeV, see Eq. (15). This suggests 
that ^H is produced in an excited state, while decaying from its ground state. We recall 
that Pauli spin is conserved in capture at rest. For K^^^^ + ^Li -^ ^H + tt"*" production, 
since ^Li is very well approximated (about 98%) by a L == 0, S = 1 configuration [26], 
^H is dominantly produced in its l"*" first excited state, decaying then by a fast magnetic 
dipole transition to the 0"*" ground state from which the mesonic weak decay occurs. 
In this situation, the pion kinetic energies T{n^) and T(7r^), directly measured by the 
FINUDA spectrometer, should reflect this systematic difference between production and 
decay. The mass of the 0+ ground state should be calculated from the decay reaction 
only, giving a mean value M(^Hg.s.) — (5800.9 ± 1.2) MeV, corresponding to a binding 



^ We thank Dr. D.J. Millener for alerting us to this estimate. A somewhat higher value, B^'^(^H) = 
4.60 ± 0.24 MeV, is obtained if the preliminary value SACXHe) = 5.68 ± 0.03(stat) ± 0.22(syst) MeV 
from the (e, e'ftr+) reaction in the JLab EOl-Oll experiment is used [27]. 
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energy of 4.5 ± 1.2 MeV with respect to (A + ^H) and of 0.8 ± 1.2 McV with respect 
to (^H + 2n). For the 1+ excited state it should be possible to evaluate a mean mass 
M(^H*) = (5801.9 ± 1.0) MeV. Although nominally unstable by 0.2 ± 1.0 MeV, the low 
Q value for two-neutron emission plus the associated AS — 1 spin flip required in the 
decay to ^H(0+ ) + 2n are likely to make the aH(1j|^j, -^ 0+^ ) Ml 7-ray transition 
competitive with the strong decay of the 1+^^, level. 

It is worth noting that the uncertainty placed on the excitation energy A_E(0+ — > 1+), 
identifying this AE with AM(^H), is considerably smaller than the uncertainty of each 
one of the 0+ and the 1+ levels because AM has been determined directly from the sum 
of kinetic energies r(7r+)+T(7r~) and its associated uncertainty. The value determined in 
the present experiment for AE is smaller by 2a than the value Aii^akaishi(0~'' — >■ 1^) = 2.4 
MeV predicted in Ref. [13]. This is in line with the conclusion drawn from the absolute 
energy location of the 0+ g.s., casting doubts on the applicability of the model developed 
by Akaishi et al. [13]. The value AEi^il : 0+ -^ 1+) == 0.98±0.74 MeV is in good accord 
with A£;(^H : 0+ -^ 1+) = 1.04±0.03 McV [2], consistently with a weak-coupling picture 
for the two 'halo' p-shcU neutrons in ^H outside the s-shell cluster of ^H. 

It is also worth noting that the width of the selected T{tt'^) + T(7r~) interval (2 McV) 
allows to include both production and decay pions within the experimental resolution 
of the mass determination at a la level, thus extending the validity of the working 
assumption on which the analysis method was based, namely that the masses of the 
produced and decaying hypernucleus are equal. A variation of the binding energy from 
the value of 5 MeV used to fix the selection on T(7r+) + T{'k^) to 4 MeV, the average 
binding energy of the 0^ and l'^ levels, produces a completely negligible variation of 
the accepted fraction of events due to the selection criteria with respect to the errors. 
Finally, the absence of systematics arising from the particular choice of both width and 
position of the selected T{'k^) -\- T{tt~) interval indicates that the difference between 
the ^H formation and decay masses is not influenced by the cut itself. It is important, 
however, to realize that the above deductions on the ^H excitation spectrum rely on very 
scarce statistics and, therefore, have to be considered as indication, even if quite solid. 

Before closing we wish to discuss briefly another scenario for the excitation spectrum 
of ^H motivated by the somewhat large spread among the three ^H candidate events. 
Apart from the 0+ g.s. and 1+ spin-flip excited state as in ^H, a 2+ excited state as for 
the p-shell dineutron system in ^He (1.80 MeV) is expected at about 2 MeV excitation in 
^H. Furthermore, a triplet of spin-flip excitations l+,2+,3+, built on the 2+ dineutron 
excitation is expected 1 MeV higher, at about 3 MeV excitation in ^H. It is then not 
unreasonable to assign event 3 in Table 1 to formation and decay of ^H involving the 
0+ g.s. and its 1+ spin-flip excited state, as considered above, whereas the other two 
events which are relatively close to each other correspond to formation of one of the 
l+,2"'",3"'" levels and to decay from the 2+ dineutron excitation. This scenario generates 
an additional excitation scale to confront the « 3 MeV separation between the first two 
events of Table 1 and the third one. This results in the assignment of ^H levels listed in 
Table 2. 

The table exhibits that the g.s. of ^H is bound in this scenario much stronger than the 
SM estimate outlined above, and in fact it is even more bound than predicted by Akaishi 
[13], although the excitation energy of the spin- flip 1+ level appears considerably smaller 
than in his prediction. The excitation energy of the 2+ level comes out about 3.4 MeV, 
too high with respect to the simple SM consideration. We conclude that this scenario is 
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Table 2 

Masses and Bj^ values (in McV) of ^H levels assuming that event 3 in Table 1 corresponds to the lowest 

two levels and events 1 and 2 correspond to higher levels. 



6H 0+ 1+ 2+ (l+,2+,3+) 

M 5798.66±0.84 5799.97±0.96 5802.07±0.59 5802.89±0.68 
Ba 6.78±0.84 5.47±0.96 3.37±0.59 2.55±0.68 

unlikely, but this conclusion does not derive from any model-independent experimental 
observation. Future experiments will tell. 

8. Conclusions 

We have reported the first observation of the hyper superheavy hydrogen ^H, based on 
detecting 3 candidate events that cannot be attributed to pure instrumental or physical 
backgrounds. The resulting binding energy of ^H, B\ = 4.0 ± 1.1 MeV, agrees with 
simple shell-model estimates initiated by Dalitz and Levi Sctti [1], but disagrees with the 
prediction made by Akaishi [13] based on a strongly attractive ANN interaction within 
a coherent AA^ — YiN mixing model. It was suggested that the excitation energy of the 
1+ spin-flip state with respect to the 0+ g.s. be identified with the systematic difference 
AM — 0.98 ± 0.74 MeV between values of ^H mass derived separately from production 
and from decay. This value is consistent with the 1.04 MeV for the analogous spin- flip 
excitation in ^H, confirming again the applicability of the shell-model estimates. An 
experiment to produce ^H via the {tt^,K^) reaction on ^Li at 1.2 GeV/c was recently 
approved at J-PARC [28] and should run soon. The expected energy resolution is 2.5 
MeV FWHM, and the expected statistics about 1-2 orders of magnitude higher than 
previous KEK experiments. 
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